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1 Introduction 

Theoretical predictions for non-leptonic decays are obtained by introducing 
an effective low-energy Hamiltonian expressed in terms of local operators 
and of the corresponding Wilson coefficients. The latter can be computed in 
perturbation theory, whereas the matrix elements of the operators have to 
be evaluated within some non-perturbative approach. For kaon decays, the 
Wilson coefficients are known at the next-to-leading order accuracy [[jj-[f§ 
and the main uncertainties come from the calculation of the matrix elements. 
In this talk we review the present status of lattice computations of matrix 
elements which are relevant in kaon decays. We focus, in particular, on 
those which enter the calculation of e'/e and make also a comparison with 
other methods, namely the Chiral Quark Model (xQM) and the large N 
expansion fl9|— 1|14| . For a more general discussion of the theory of CP violation 



in kaon decays see |15[ and references therein. 

Some general remarks are necessary before entering a more detailed dis- 
cussion. Given the large numerical cancellations which may occur in the 
theoretical expression of e'/e, a solid prediction should avoid the "Harlequin 
procedure" . This procedure consists in patching together B 6 from the xQM, 
B 8 from the 1/N expansion, m^ 5 from the lattice, etc., or any other combi- 
nation/average of different methods. All these quantities are indeed strongly 
correlated (for example B 6 and B 8 in the 1/N expansion or B parameters 
and quark masses in the lattice approach [16|]) and should be consistently 
computed within each given theoretical framework. Unfortunately, none of 
the actual non-perturbative methods is in the position to avoid completely 
the Harlequin procedure, not even for the most important input parameters 
only. The second important issue is the consistency of the renormalization 
procedure adopted in the perturbative calculation of the Wilson coefficients 
and in the non-perturbative computation of the operator matrix elements. 
This problem is particularly serious for the xQM and the 1/N expansion, and 
will be discussed when comparing the lattice approach to these methods. We 
will address, in particular, the problem of the quadratic divergences appear- 
ing in the 1/N expansion. This is an important issue, since the authors of 
ref. [0, [L4| find that these divergences provide the enhancement necessary 
to explain the large values of ReA Q and of e'/e. 



Schematically, e' can be cast in the form 

, exp(m/4) uj 



x 



oj~ 1 ImA 2 - (1 - Vt IB )ImA (1) 



V2 ReA 
where u = ReA 2 / 'ReA and ReA are taken from experiments, and Qjb is a 



correcting factor, estimated in refs. [17|— [19|, due to isospin-breaking effects. 



Using the operator product expansion, the K — > im amplitudes ImA 2 and 
ImAo are computed from the matrix elements of the effective Hamiltonian, 
expressed in terms of Wilson coefficients and renormalized operators 

(7in\H AS=1 \K°) I=0 , 2 = J2(^\Q^)\K°) 1=0,2 Q(fi) (2) 

i 

where the sum is over a complete set of operators, which depend on the 
renormalization scale /i. Wilson coefficients and matrix elements of the op- 
erators Qi(fi), appearing in the effective Hamiltonian, separately depend on 
the choice of the renormalization scale and scheme. This dependence cancels 
in physical quantities, such as ImA 2 and ImAo, up to higher-order cor- 
rections in the perturbative expansion of the Wilson coefficients. For this 
crucial cancellation to take place, the non-perturbative method used to com- 
pute hadronic matrix elements must allow a definition of the renormalized 
operators consistent with the scheme used in the calculation of the Wilson 
coefficients. 

So far, lattice QCD is the only non-perturbative approach in which both 
the scale and scheme dependence can be consistently accounted for, using 
either lattice perturbation theory or non-perturbative renormalization tech- 
niques |^0|, |21|. This is the main reason why the authors of refs. [||, ^2], ^3], ^4 



have followed this approach over the years. 

There is a general consensus |T3 that the largest contributions are those 



coming from Qq and Qs (for ImA 2 ), with opposite sign, and sizeable con- 
tributions may come from Q3, Q4 and Qg (for ImA 2 ) in the presence of 
large cancellations between Q 6 and Q&, i.e. when the prediction for e'/e is 
of O(10~ A ) Q. For this reason the following discussion, and the comparison 
with other calculations, will be focused on the determination, and errors, of 
the matrix elements of the two most important operators. 



lr rhe operator (33,4,5,6 only contribute to the 1 = amplitudes. 



2 Matrix elements from lattice QCD 

The evaluation of physical K — > tttt matrix elements on the lattice relies 
on the use of Chiral Perturbation Theory (%PT): so far only (ir\Qi(fj,)\K) 
and (tt(p = 0)ir(q = 0)\Q i (fj,)\K)i =2 (with the two pions at rest) have been 
computed for a variety of operators. The physical matrix elements are then 
obtained by using %PT at the lowest order. This is a consequence of the dif- 
ficulties in extracting physical multi-particle amplitudes in Euclidean space- 
time P5| . Proposals to overcome this problem have been presented, at the 
price of introducing some model dependence in the lattice results lp6fl . The 
use of x?T implies that large systematic errors may occur in the presence of 
large corrections from higher-order terms in the chiral expansion and/or from 
FSI. This problem is common to all approaches: if large higher-order terms 
in the chiral expansion are indeed present and important, any method aiming 
to have these systematic errors under control must be able to reproduce the 
FSI phases of the physical amplitudes. The approaches of ref. R ITDJ and 



[|T|, [13|, [14], J27J, however, give FSI smaller than their physical values. 



1 = 2 matrix element of Q%. There exists a large set of quenched calcula- 
tions of {Q&)2 performed with different formulations of the lattice fermion ac- 
tions (Staggered, Wilson, tree-level improved, tadpole improved) and renor- 
malization techniques (perturbative, boosted perturbative, non-perturbative) 
at several values of the inverse lattice spacing a -1 =2-^3 GeV JL6], ^1| ^8], |9[ 
30fl . All these calculations, usually expressed in terms of B% , give consis- 



tent results within 20% of uncertainty. Among the results, in the numerical 
estimates presented in sect. [|, we have taken the central value from the re- 



cent calculation of ref. [IB|, where the matrix elements {Q&)2 and ((57)2 have 
been computed without any reference to the quark masses, and inflated the 
errors to account from the uncertainty due to the quenched approximation 
(unquenched results are expected very soon) and the lack of extrapolation to 
zero lattice spacing. For a discussion on the role of quark masses see below 
at the end of this section. 

Matrix element of Qq. For (Qq)q from the lattice, the situation appears 
worse today than a few years ago when the calculations of refs. [§, |22|, p3| 
were performed: 

i) until 1997, the only existing lattice result, obtained with staggered 



fermions (SF) without NLO lattice perturbative corrections, was B$ = 1.0 ± 
0.2 |3lH - This is the value used in previous analyses || ^2], p3| ; 

ii) with SF even more accurate results have been quoted recently, namely 
B 6 = 0.67 ± 0.04 ± 0.05 (quenched) and B & = 0.76 ± 0.03 ± 0.05 (with 

n f = 2) H; 

iii) 0(a s ) corrections, necessary to match lattice operators to continuum 
ones at the NLO, are so huge for Q$ in the case of SF (in the neighbourhood 
of —100% [P^D as to make all the above results unreliable. Note, however, 
that the corrections tend to diminish the value of (Qe)o; 

iv) the latest lattice results for this matrix element, computed with 
domain-wall fermions [|3|] from (tt\Qq\K), are absolutely surprising: (Qq)o 



has a sign opposite to what expected in the VSA, and to what is found with 
the xQM and the 1/N expansion. Moreover, the absolute value is so large 
as to give e'/e ~ —120 x 10 -4 . Were this confirmed, even the conservative 
statement by Andrzej Buras fL5| , namely ... that certain features present 
in the Standard Model are confirmed by the experimental results. Indeed the 
sign and the order of magnitude of e'/e predicted by the SM turn out to agree 
with the data... would result too optimistic. In order to reproduce the ex- 
perimental number, e'/e ~ 20 x 10~ 4 , not only new physics is required, but a 
large cancellation should also occur between the Standard Model and the new 
physics contributions. Since this result has been obtained with domain-wall 
fermions, a lattice formulation for which numerical studies started very re- 
cently, and no details on the renormalization and subtraction procedure have 
been given, we consider premature to use the value of the matrix element 
of ref. |34] in phenomenological analyses. Hopefully, new lattice calculations 



will clarify this fundamental issue. 

B parameters and quark masses. Following the common lore, matrix el- 
ements of weak four-fermion operators are given in terms of the so-called 
^-parameters which measure the deviation of their values from those ob- 
tained in the Vacuum Saturation Approximation (VSA). A classical exam- 
ple is provided by the matrix element of the AS = 2 left-left operator 
Q AS=2 = S7 M (1 — 7 5 )<i S7 M (1 — 75)6? relevant to the prediction of the CP- 
violation parameter e 

(K°\Q* s =z\K°) = \M\f K B K . (3) 



VSA values and 5-parameters are also used for matrix elements of operators 
entering the expression of e'/e, in particular Q 6 = 5^7^(1 — 75)^/3 J2q 9/37^(1 + 
7 5 )g Q and Q 8 = 3/2^7^(1 - 75)^ Y, q e q q^{l + ls)<la 



(7nr\QM\K ) 
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m a (/x) +m d (ji) 



(nn\Q 8 (Li)\K ) I=2 = yfif* 
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(f K -U)BM 
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- [Ml - Ml 



m s (fi) +m d (jj,)J 
Bf l2 \») • 



(4) 



Since in the VSA and in the 1/N expansion the expression of the matrix 
elements is quadratic in m s + m^, predictions for the physical amplitudes 
are heavily affected by the specific value which we assume for this quantity. 
Contrary to fx, M^, quark masses are not directly measured by experiments 
and the present accuracy in their determination is still rather poor [[!5|, |36| 



Therefore, the "conventional" parametrization induces a large systematic un- 
certainty in the prediction of the physical amplitudes of (Qe)i=o and (Qs)i=2 
(and of any other left-right operator). Moreover, whereas for Q AS=2 we in- 
troduce Bk as an alias of the matrix element, by using (W we replace each 
of the matrix elements with 2 unknown quantities, i.e. the 5-parameter 
and m s + rrid- Finally, in many phenomenological analyses, the values of 
the 5-parameters of (Qe)i=o and (Qs)i=2 and of the quark masses are taken 
by independent lattice calculations, thus increasing the spread of the theo- 
retical predictions. All this can be avoided in the lattice approach, where 
matrix elements can be computed from first principles. In ref. |16J a new 
parametrisation of the matrix elements in terms of well known experimen- 
tal quantities, without any reference the strange (down) quark mass, has 
been introduced. This results in a determination of physical amplitudes with 
smaller systematic errors. The interested reader can refer to [HJ for details. 

Before ending this discussion, we wish to illustrate the correlation existing 
between the B parameters and the quark masses in lattice calculations. On 
the lattice, quark masses are often extracted from the matrix elements of 
the (renormalized) axial current {A^) and pseudoscalar density (P(/x)) (for 



simplicity we assume degenerate quark masses) 

m(/i) -2(a|P(^)|/5)' (5) 

where a and (3 are physical states (typically a is the vacuum state and (3 
the one-pseudoscalar meson state) and m(n) and P(/x) are renormalized in 
the same scheme. On the other hand, the B parameters of Qq and Q$ are 
obtained (schematically) from the ratio of the following matrix elements, 
evaluated using suitable ratios of correlation functions 0: 

B * Moc (n\PM\o)(o\PM\K)> (6) 

where P^ and Pjc are the pseudoscalar densities with the flavour content 
of the pion or kaon, respectively. Eqs. @ and (|6]) demonstrate the strong 
correlation existing between B parameters and quark masses: large values of 
the matrix elements of P{^) correspond, at the same time, to small values 
of ra(/x) and i?6,8 (/•*)• Physical amplitudes, instead, behave as 

(Q 6 , 8 ) = const, x ^M (7 ) 

where "const." is a constant which may be expressed in terms of measurable 
quantities (specifically M K and fx) only. From eqs. (§) and (^|), we recognize 
that the dependence on (P(/x)) cancels in the ratio 5 6j8 /m(/i) 2 , appearing in 
the physical matrix elements. 

Previous lattice studies preferred to work with B parameters because 
these are dimensionless quantities, not affected by the uncertainty due to the 
calibration of the lattice spacing. This method can still be used, provided 
that quark masses and the B parameters from the same simulation are pre- 
sented together (alternatively one can give directly the ratio B 6<8 /m(/i) 2 ). 
In ref. |[16|| , two possible definitions of dimensionless li B parameters", which 



can be directly related to physical matrix elements without using the quark 
masses have been proposed. 

The strange quark mass. Although in lattice calculations of matrix ele- 
ments any reference to quark masses can be avoided, these are fundamental 

2 See for example ref. pl[. For simplicity the superscript (3/2) in _Bg is omitted. 



parameters of the Standard Model and are used in the large N expansion. 



Here we would like to add only a few remarks to ref. |36| , where this subject 
has been reviewed. 

The first observation is the following. Lattice calculations that use non- 
perturbative renormalization methods obtain the quark masses without er- 
rors coming from the truncation of the perturbative series (typically in the 
RI-MOM or the Schrodinger renormalization schemes; for a complete set of 
references see ref. 113611). The conversion of these results to the "standard" 



MS scheme can be done at the N 3 LO. Differences between NLO, N 2 LO and 
N 3 LO are important, ~ 6 -f- 10 MeV for m^ s , as demonstrated by the follow- 
ing example 

NLO N 2 LO N 3 LO 

mp(2GeV) = {5.2(5); 4.9(5); 4.8(5)} MeV, 

mp(2GeV) = {120(9); 114(9); 111(9)} MeV, (8) 

taken from ref. [[J7[]. Therefore, when confronting results from different cal- 
culations it is necessary to specify the order at which the results have been 



obtained. In table 1 of ref. |5B|], results obtained with perturbation theory at 
NLO or with non-perturbative methods at the N 3 LO are directly compared. 
This, for the reasons discussed above, is misleading. We also note that in 
most of the phenomenological applications, for example with QCD sum rules, 
the theoretical expressions are only known at the NLO and, for consistency, 
quark masses at the same level of accuracy should be used. 

By comparing the results obtained with the non-perturbatively improved 
action at /3 = 6.2 (corresponding to a" 1 ~ 2.6 GeV, which is the value used by 
the APE Collaboration) with those extrapolated to the continuum (table 2 
of ref. ||38||), one finds the discretization errors at this value of the lattice 
spacing and with this action to be 3 -j- 4%. This is much smaller than the 
15% quoted in BH . 



The large reduction of the value of the masses in the unquenched case, 
found by the CP-PACS Collaboration, is not confirmed by other lattice cal- 



culations by the MILC [[39| and APE Collaborations HDJ and it is at variance 
with the bounds of ref. ||l| . We think that further investigation is required 
on this important point. 
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3 Renormalization group invariant operators 

Wilson coefficients and renormalized operators are usually defined in a given 
scheme (HV, NDR, RI), at a fixed renormalization scale /i, and depend on 
the renormalization scheme and scale. This is a source of confusion in the 
literature. Quite often, for example, one finds comparisons of B parameters 
computed in different schemes. Incidentally, we note that the NDR scheme 
used in the lattice calculation of ref. |29| differs from the standard NDR 



scheme of refs. 0-0; on the other hand, the HV scheme of ref. || is not 
the same as the HV scheme of ref. @. In some cases, the differences between 
different schemes may be numerically large, e.g. B\ ~ 1.3 B% 

at fi ~ 2 GeV. To avoid all these problems, it is convenient to introduce a 
Renormalization Group Invariant (RGI) definition of Wilson coefficients and 
composite operators which generalises what is usually done for B K using the 
RGI 5-parameter Bk- The idea is very simple. Physical amplitudes can be 
written as 

(F\H\I) = {F\Q(»)\I).C(»), (9) 

where Q(fi) = (Qi (/•*), Q 2(^)1 ■ ■ ■ , Qn(i^)) is the operator basis and C(n) the 
corresponding Wilson coefficients, represented as a column vector. C(n) is 
expressed in terms of its counter-part, computed at a large scale M, through 
the renormalization-group evolution matrix W\p, M] 

d(ji) = W\p,M\C(M). (10) 

The initial conditions for the evolution equations, C(M), are obtained by 
perturbative matching of the full theory, which includes propagating heavy- 
vector bosons (W and Z°), the top quark, SUSY particles, etc., to the ef- 
fective theory where the W, Z°, the top quark and all the heavy particles 
have been integrated out. In general, C(M) depends on the scheme used to 
define the renormalized operators. It is possible to show that W\ja, M] can 
be written in the form 

W[fi,M] = M[^]U[^M}M- l [M] , (11) 

with 

\a (M)V^ m . * a(u) * 

U\fi,M\= ^TT > %H1 + T^«], (12) 



(n\ r p A 

where 7q is the leading order anomalous dimension matrix and J[A(/i)] can 
be obtained by solving the Renormalization Group Equations (RGE) at the 
NLO. By defining 

tS-fy] = MM [a a (n)r$ )T/2p0 , (13) 

we get 

W[fj,,M] =w~ 1 [/j]w[M} . (14) 

The effective Hamiltonian @ can then be written as 

ft = Q(fx)-C(fx)=Q^)W[^M]C(M) 

= Qi^w- 1 ^} ■ w[M}C(M) = Q RGI ■ C RGI , (15) 

with 

C RGI = w[M)C(M) , Q RGI = Q{n)-w~ l [n\ . (16) 

C and Q^ 07 are scheme and scale independent at the order at which the 
Wilson coefficients have been computed. 

4 Comparison with other methods 

In this section, we briefly discuss the relevant aspects which distinguish the 
lattice approach from others which have been used in the literature to predict 

e'/e. 

The original approach of the Munich group was to extract the values of 
the relevant matrix elements from experimental measurements [|], |B[|. This 
method guarantees the consistency of the operator matrix elements with the 
corresponding Wilson coefficients. Unfortunately, with the Munich method 
it is impossible to get the two most important contributions, namely those 
corresponding to (Qe)o an d {Qs}2- For this reason, "guided by the results 
presented above and biased to some extent by the results from the large- 



iV approach and lattice calculations", the authors of ref. ]42| have taken 
B 6 = 1.0 ± 0.3 and B^ 3/2) = 0.8 ± 0.2, and mf^ = 110 ± 20 MeV at n = 1.3 
GeV. These values, if assumed to hold in the HV regularization, are very 
close to those used in ref. [24| . They do not come however from a calculation 



consistently made within a given theoretical approach (large N expansion, 
xQM or lattice for example). 

The 1/N expansion and the xQM are effective theories. To be specific, 
in the framework of the 1/N expansion the starting point is given by the 
chiral Lagrangian for pseudoscalar mesons expanded in powers of masses 
and momenta. At the leading order in 1/N, local four-fermion operators can 
be written in terms of products of currents and densities, which are expressed 
in terms of the fields and coupling of the effective theory. In higher orders, a 
(hard) cutoff, A c , must be introduced to compute the relevant loop diagrams. 
The cutoff is usually identified with the scale at which the short-distance 
Wilson coefficients must be evaluated. 

Divergences appearing in factorizable contributions can be reabsorbed 
in the renormalized coupling of the effective Lagrangian and in the quark 
masses, non-factorizable corrections constitute the part which should be 
matched to the short distance coefficients. By using the intermediate colour- 
singlet boson method, the authors of refs. |Tl|]-|l4|],[[|3l claim to be able to 



perform a consistent matching, including the finite terms, of the matrix el- 
ements of the operators in the effective theory to the corresponding Wilson 
coefficients. It is precisely this point which, in our opinion, has never been 
demonstrated in a convincing way. If the matching is "consistent" , then it 
should be possible to show analytically that the cutoff dependence of the 
matrix elements computed in the 1/N expansion cancels that of the Wilson 
coefficients, at least at the order in 1/N at which they are working. Moreover, 
if really finite terms are under control, it should be possible to tell whether 
the coefficients should be taken in HV , NDR or any other renormalization 
scheme. 

The fact that in higher orders even quadratic divergences appear, with 
the result that the logarithmic divergences depend now on the regularization, 
makes the matching even more problematic. Theoretically, we cannot imag- 
ine any mechanism to cancel the cutoff dependence of the physical amplitude 
in the presence of quadratic divergences, which should, in our opinion, dis- 
appear in any reasonable version of the effective theory. It is also important 
to show (and to our knowledge it has never been done) that the numerical 
results for the matrix elements are stable with respect to the choice of the 
ultraviolet cutoff. This would also clarify the issue of the routing of the mo- 



ld 



menta in divergent integrals. For example, the matrix elements in the meson 
theory could be computed in some lattice regularization. 



5 Numerical results 

As discussed above, all the methods used in the calculation of e'/e are not 
completely satisfactory and in general suffer from large theoretical uncertain- 
ties. 

The lattice approach can, in principle, compute the relevant matrix el- 
ements without any model assumption (at least at the lowest order in the 
chiral expansion), and with operators consistently defined to match the Wil- 
son coefficients of the effective Hamiltonian. In spite of these advantages 
the lattice results for (Qq)q are inconclusive, as discussed before. Regarding 



the surprising result of ref . [g4| , we think that further scrutiny and confirma- 
tion from other calculations are needed before using it in a phenomenological 
analysis. 



In the absence of any definite result for (Qe)o from the lattice, ref. p4 
assumed for this matrix element the value 

(Q 6 ) = (TrvrlQf V \K°) I=0 = -0.4 ± 0.4 GeV 3 , (17) 

and 

(Q 5 >o = l/3(Q 6 >o, (18) 

at a scale fi = 2 GeV. The value of the matrix element in eq. (|T7|) corresponds 
to B % = l.Oil.O for a "conventional" mass fixed to mf IS + mf s = 130 MeV. 

For (Qt ) 8)2, the values of ref. [[U| (obtained with an improved action 



using non-perturbatively renormalized operators at /i = 2 GeV) have been 
used, namely 

(Q 7 ) 2 = (nn\Q? v \K°) I=2 = 0.18 ± 0.06 GeV 3 , (19) 

(Q 8 ) 2 = (nn\Qf v \K°) I=2 = 0.62 ± 0.12 GeV 3 , (20) 

where the superscript HV denotes the t'Hooft-Veltman renormalization scheme. 
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By varying the input parameters as described in ref. [El] and by weighting 



the Monte Carlo events with the experimental constraints, the prediction for 
e'/e is 

e'/e = (3.6ljj;£ ± 0.5) x 10" 4 , (21) 

where the third error on e'/e is an estimate of the residual scheme depen- 
dence due to unknown higher-order corrections in the perturbative expan- 
sion []. Given the large theoretical uncertainties, and taking into account 
some differences in the calculation of this quantity (choice of the renormal- 
ization scale, values of several B parameters, etc.), the result in eq. (|21|) is 
in substantial agreement, though slightly lower, with the recently upgraded 
evaluation of ref. g|: e'/e = (7.7±|:g) x 10" 4 and e'/e = (5.2±£?) x 10" 4 
in NDR and in HV respectively. It is also very close to previous estimates 



of the Rome [E_S, Ej| and Munich group |4], [J . This agreement it is not sur- 
prising since the two groups use very similar inputs for the matrix elements 
and the experimental parameters have only slightly changed in the last few 
years. The crucial question, namely a quantitative determination of {Qe)o, 
remains unfortunately still unsolved. 

All the above results are, however, much lower than the recent mea- 
surements of KTeV, Re(e'/e) = (28.0 ± 4.1) x 10" 4 , of NA48, Re(e'/e) = 
(18.5 ± 7.3) x 10~ 4 , or than the present world average Re(e'/e)wA = (21.2 ± 
4.6) x 10~ 4 , determined by the results of refs. p4||-|f46 



By scanning various input parameters (in the conventional approach B e , 
Bg ' , a s (M z ), Im\ t etc.) and in particular by choosing them close to their 
extreme values it is possible to obtain e'/e ~ 20 x 10 -4 . This also gives the 
impression of a better agreement (lesser disagreement) between the theoret- 
ical predictions and the data. For example, in ref. |24| the scanning gives 
— 11 x 10~ 4 < e'/e < 27 x 10~ 4 . In spite of the fact that the experimental 
world average is compatible with the "scanned" range above, a conspiracy of 
several inputs in the same direction is necessary in order to get a large value 
of e'/e. For central values of the parameters, the predictions are, in general, 
much lower than the experimental results. For this reason, barring the pos- 
sibility of new physics effects [ |48|| , we believe that an important message is 

3 The value in eq. O is slightly different from that presented at the the Conference and 
quoted by ref.Jl^]. The reason is that the final analysis of ref. Q found for (Qs)i=2 a 
value larger by about 15% than the preliminary one. 
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arriving from the experimental results: 

penguin contractions (or eye-diagrams, not to be confused with penguin 
operators), which are usually neglected within factorization, give contribu- 
tions which makes the matrix elements definitely larger than their factorized 
values. 

This implies that the "effective" B parameters of the relevant operators, 
specifically those relative to the matrix elements of Qi and Q2 for Re(Ao) and 
of Qq for e'/e are much larger than 1. This interpretation would provide a 
unique dynamical mechanism to explain both the AI = 1/2 rule and a large 



value of e'/e [49]. Large contributions from penguin contractions are actually 
found by calculations performed in the framework of the Chiral Quark Model 
(xQM) |9|1-[ r0|1 or the 1/N expansion [11, [H| 14), |27]. It is very important 



that these indications find quantitative confirmation in other approaches, 
for example in lattice QCD calculations. Note that naive explanations of 
the large value of e'/e, such as a very low value of m* /s , would leave the 
AI =1/2 rule unexplained. 

Finally, one may try to quantify the amount of enhancement required for 
the matrix element of Qq in order to explain the experimental value of e'/e. 
A fit of (Q 6 ) to Re(e'/e) WA gives (Q 6 ) = — l-2±g;lf ± 0.15 GeV 3 , about 2.5 
times larger than the central value used in our analysis. 
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